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We study the entangling power of a nanoelectromechanical system (NEMS) simultaneously interacting with 
two separately trapped ions. To highlight this entangling capability, we consider a special regime where the ion- 
ion coupling does not generate entanglement in the system, and any resulting entanglement will be the result 
of the NEMS acting as an entangling device. We study the dynamical behavior of the bipartite NEMS-induced 
ion-ion entanglement as well as the tripartite entanglement of the whole system (ions-l-NEMS). We found some 
quite remarkable phenomena in this hybrid system. For instance, the two trapped ions initially uncorrelated 
and prepared in coherent states can become entangled by interacting with a nanoelectromechanical resonator 
(also prepared in a coherent state) as soon as the ion-NEMS coupling achieve a certain value, and this can be 
controlled by external voltage gate on the NEMS device. By considering the NEMS in an initial thermal state, we 
numerically show that there is not a temperature threshold above which bipartite ion-ion entanglement ceases. 
A distinct effect occurs when the NEMS interacts with a thermal reservoir, above a certain value of temperature, 
the NEMS induction of ion-ion entanglement ceases. We also show that tripartite entanglement presents a more 
pronounced robustness against the destructive effects of dissipation when compared to the bipartite content. 

PACS numbers: 03.67.Bg,85.85.+j 



I. INTRODUCTION 

During the last years, theoretical and instrumental develop- 
ments in physics have led to an unprecedented level of con- 
trol over atomic and nanoscale systems. Striking examples 
include the manipulation of single quantum systems such as 
laser cooled trapped ions |T|, and the advances on fabrica- 
tion, characterization and application of nano- mesoscale de- 
vices such as electromechanical resonators ||2l- Isolated single 
quantum systems are the natural choice when quantum logic 
is involved |3| and the nanoscale devices are recognizedly 
known as very sensitive weak signal detectors [4-6J. The 
emerging field of quantum technology takes advantage of the 
strong and different features of these two kind of systems by 
combining them together in a single setup |7-9|. 

The goal of quantum technology is to explore legitimate 
quantum resources such as entanglement to perform tasks 
with no classical analogue and to enforce it on nano- meso- 
scopic systems which can be integrated, for instance, on a 
chip for massive use. In order to achieve that, one very 
attractive route consists of cooling the NEMS down to its 
ground state |10| and then enforce on it a legitimate quan- 
tum behavior by coupling it to well controlled small quantum 
systems. For example, the coupling of qubits to these me- 
chanical nanoresonators allows one to generate superpositions 
of macroscopically distinguishable states in the motion 111]. 
These states may be useful to study the so called quantum-to- 
classical transition and the phenomenon of decoherence hap- 
pening in nano- and mesosystems. This is a natural follow up 
to the studies started in the scope of cavity quantum electro- 
dynamics |12|. On the other hand, potential applications of 
NEMSs include several uses as a sensitive detector Detec- 
tion of single spins |4| and spin-spin interactions |5 | as well 
as single molecule mass spectrometry lfT3l offer new possibil- 
ities in chemical imaging and characterization. For instance. 



nitrogen-vacancy centers in a diamond tip coupled to NEMS 
resonators operating at room temperature 1 14] offers a new 
route to detect a chemical element that carry a nuclear mag- 
netic moment. This device is sensitive enough to determine 
the element identity and arrangement in a complex molecule. 
In the case of mass spectroscopy, NEMS resonators are be- 
lieved to enable in the future a detection device operating with 
resolution better than 1 Dalton (the hydrogen atom mass) 1 15 1. 
Nowadays, the state-of-the-art is that single gold atoms can be 
measured I.15J. 

Trapped ions are well studied quantum systems where 
experimental control of coherent quantum dynamics, state 
preparation, and measurement are achieved with great pre- 
cision and good fidelity |[T]. As a matter of fact, trapped 
ions have been demonstrated to achieve an error probability 
per randomized single-qubit gate below the threshold esti- 
mate commonly considered sufficient for fault-tolerant quan- 
tum computing |16|. This achievement of low single -qubit- 
gate errors is an essential step toward scalable quantum com- 
puters and simulators. It is then clear that a hybrid system 
composed of trapped ions and NEMS resonators represent an 
interesting platform to study the interplay of quantum me- 
chanical effects in atomic and nanoscopic systems. Follow- 
ing this idea, one proposal L7J suggest the use of a trapped 
ion as a probe or device controller acting on the NEMS. In 
another proposal fE], trapped ions are applied to monitor and 
manipulate the number state of a NEMS resonator such that it 
enables a statistical inference of the mean phonon number of 
the oscillator. 

Here, we take a step forward and propose the use of 
a NEMS resonator to generate entanglement between two 
trapped ions. This is a conceptually important step since it 
goes the other way around of what is normally performed. In- 
stead of using the microsystem to cause quantum behavior of 
the mesoscopic system, we use the latter to enforce quantum 
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correlations on the former Since a classical resonator is un- 
able to generate entanglement among quantum systems cou- 
pled to it, this kind of study may work as a kind of classicality 
test in what concerns entanglement generation. We also study 
the qualitative behavior of tripartite entanglement as given by 
residual tripartite entanglement ifTTl and contrast its dynam- 
ical behavior to the ones of the bipartite entanglements li 1 8i . 
Moreover, we investigate the regime where NEMS is initially 
in a thermal state, in such case one could expect that bipar- 
tite ion-ion entanglement would die for high temperature, but 
surprisingly there is not a temperature threshold above which 
they do not get entangled. Finally, we present results that 
show that in this system tripartite entanglement is more robust 
against energy relaxation than the bipartite entanglement. 

The paper is organized as follows. Section|ll]is dedicated to 
the presentation of the physical system studied here and Sec- 
tion III presents the theoretical tools necessary to understand 



the results presented in Section|lV] Finally, Section pTV D| con- 
tains a summary of our results and conclusions. A appendix 
shows details about lengthy analytical expressions. 



n. PHYSICAL SYSTEM AND MODEL 

In this work, we consider a system consisting of two sep- 
arate trapped ions capacitively coupled to a small doubly 
clamped nano beam (a NEMS resonator). The system is de- 
picted in figure [T] The ions-NEMS coupling results from the 
application of an external bias voltage at an electrode on the 
NEMS (the bias gate is not shown in the picture). 




Figure 1 : Sketch of the physical setup considered in this paper. Two 
trapped ions (traps not shown) interact with the flexural motion of 
a doubly clamped nano beam which is charged by a bias gate not 
shown in the picture. 

The interaction energy between the motional degree of free- 
dom of an ion with charge +e and a NEMS with bias gate 
charge Q is given by the electrostatic coupling Q [H 
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where Xi and x refer to position operators of one ion and the 
NEMS, respectively, around their equilibrium positions. At 



equilibrium, each ion is separated from the center of mass po- 
sition of the NEMS by a distance d. Since we are dealing with 
positions around the equilibrium, the signs ± are necessary to 
distinguish one ion from the other, for example, one can asso- 
ciate + with ion 1 and — with ion 2. Typically, one finds that 
{Xi), (x) << d so that the coupling energy can be expanded 
as 
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with Q = CqVo, where Cq and Vo are the capacitance and 
gate voltage, respectively. According to|2] the coupling be- 
tween the ifh ion and the NEMS is just V — ~x^iX, where 
X — . The linear and quadratic terms involving Xi or 

X alone have been absorbed into redefinitions of equilibrium 
positions and oscillators' frequencies. 

Since the ions are not neutral particles, they will also cou- 
ple by means of the same mechanism, i.e. and interaction 
term proportional to X1X2 lfT9]| . However, once the equilib- 
rium distance between them is 2d (twice the distance to the 
NEMS) and we will be considering identical traps (identical 
oscillation frequencies), only resonant terms in X1X2 will be 
relevant in this weak coupling regime (rotating wave approx- 
imation RWA) \19}. Considering V — ~x-^i^ the inter- 
action between each ion and the NEMS and a weak RWA in- 
teraction between the ions, the Hamiltonian operator for this 
tripartite mechanical oscillating system reads (ft = 1) 



H 



waJflQ + Vja}-dj — ^{ai^di + d\d2) 

J = l,2 



(3) 



where ao(aQ) and ai{a\) are the annihilation(creation) oper- 
ators for the NEMS and ion i, respectively. Also, Vi is the 
vibrational frequency of each ion (trap frequency), Vl is the 
RWA interaction strength between the ions, m; and M the 
masses of the ith ion and NEMS, respectively, and 
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X 



(4) 



is the coupling constant for interaction between the NEMS 
and the ith ion. 

It is worthwhile to notice that the vibrational motion of the 
ions are directly coupled by means of a beam-splitter-like in- 
teraction which is well known to not dynamically generate 
entanglement for direct products of thermal or coherent states 
EOl . In this situation, by putting such ions in interaction with 
the NEMS, initial creation of entanglement will be the result 
of the action of this third party acting as an entangling de- 
vice. The key element for the initial creation of ion-ion en- 
tanglement is the presence of terms of the kind dodi or d\d\ 
with strength providing an indirect squeezing interaction 
between the ions. This is an interesting problem because, be- 
ing an indirect interaction, not all values of Kj will lead to 
squeezing and entanglement in the ionic subsystem. Studying 
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the conditions for this to happen is one of the goals of this 
article. As a final remark, it is important to emphasize that 
Hamiltonian Q is quadratic in the creation and annihilation 
operators and then preserves Gaussianity. 



III. THEORETICAL TOOLS 

In this section, we present the main tools used in this work. 
Basically, we are going to discuss some techniques to treat 
continuous variable systems in Gaussian states and also mea- 
sures of entanglement suitable to study tripartite systems. 



A. Gaussian Systems 

In order to study the entanglement dynamics that results 
from the quadratic Hamiltonian (|3]l, we represent it in terms 
of coefficients for the momenta and coordinate operators 
grouped together in the vector R — {xq,pq,xi,pi,X2,P2)^ 
which contains momentum and position operators ordered in 
a convenient manner. The result is 



with 



i,j=0 



(6) 



where ( • ) denotes an expectation value. For a bipartite system 
composed of subsystems A and B [301 , an entanglement mea- 
sure called logarithmic negativity can be written as a function 
of the covariance matrix as 121] 



A\B 



= -i^ln[min(l,Ai7) 



(10) 



where 11^ are the symplectic eigenvalues fST] of the co- 
variance matrix 7^^, evaluated after partial transposition of 
system B. This partial transposition is achieved by a local 
time inversion in the oscillators pertaining to subsystem B: 
p — > —p. For example, if the system AB is composed by two 
oscillators, 7^^ = V-fV, with 



7' = Diag(l, 1,1,-1). 



(11) 



In order to study the time evolution of entanglement in this 
system, we make explicit use of the fact that (|3]l preserves the 
Gaussian character of the global state, so that we just have to 
calculate 'y{t). For a system Hamiltonian of the form (jsjl, the 
CM evolves as 
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1 1 ^ 

H = -R^nk+K =^Y. {i^^^o^J+P^^^^oPJ)+K, (5) where 



Et = exp [mt] 



(12) 



(13) 



is a symplectic matrix obeying E^JEj = J. As a result, all 
closed system dynamics is dictated by the eigenvalues of iT-L. 



B. Tripartite Entanglement 




(7) 



and, finally, K = \{uj + vi + V2) \s a. constant. The operators 
composing the vector R fulfill the usual canonical commuta- 
tion relation: Rk\ — i-ijk, where the symplectic (6 x 6) 
diagonal-block matrix J is given by: 



(8) 



In order to explore the entanglement content of the present 
model, we present now (and in next subsection) the neces- 
sary mathematical formalism of Gaussian CV entanglement 
theory. 

Choosing initial Gaussian states for the three oscillators, 
Hamiltonian (j5]l will preserve Gaussianity and we can extract 
all information about the system from the knowledge of the 
accompanying covariance matrix (CM) 7 whose elements are 
given by 
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^{RjRk + RkRj) — {Rj){Rk) 



(9) 



Our main goal in this work is to study the dynamics of en- 
tanglement in the present system, especially because it is gen- 
erated by the coupling of a nanoscale system to two atomic 
ions, forming an interesting platform for quantum technology 
experiments [Z ^. Bipartite entanglement can be quantita- 
tively addressed in a Gaussian system by using the tools devel- 
oped in the last section. We now want to study the build up of 
genuine tripartite entanglement in this hybrid system. If one 
wants to address this issue quantitatively, the way is to start 
from the well known monogamy relation for pure states fiTl 
and define a legitimate tripartite measurement called residual 
tripartite entanglement by using concurrence or logarithmic 
negativity. In the case of mixed states, this quantification has 
to be done by means of a convex-roof extension ll24l . 

For the purposes of our investigation, a semiquantitative 
approach is sufficient since we are not aiming at a particu- 
lar quantum information protocol where a given value of an 
entanglement measure has a precise meaning. Instead, we are 
interested in knowing whether or not there is legitimate tripar- 
tite entanglement in our system when some bipartite entan- 
glement has vanished and also to investigate the robustness of 
bi- or tripartite entanglement against dissipation in the NEMS. 
We can then follow the same procedure taken in |25 1 which is 
based on multimode inseparability classification proposed by 
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Giedke et al. 11261 . In our case this would correspond to check 
whether or not the negativity ( [T0| vanishes for all bipartitions, 
i.e., by assigning {0} to the NEMS and {1, 2} to the ions, we 



as well as 



want to check when N^j^. =^ for all k — 0, 1, 2 with 
* 7^ J 7^ ^ — according to L26J, the system is genuinely tri- 
partite entangled only if this is the case. Consequently, one 
can, for instance, study the resilience of tripartite entangle- 
ment against dissipation by investigating the quantity 



roi2 = min[A'"o|i2,^i|02,^^2|oi] 



From the multimode inseparability classification |26| it is then 
clear that genuine tripartite entanglement ceases to exist as 
soon as roi2 = 0. 



Al = y|(iV + l)(i,-l)^ and \2 = -\j\N 

(20) 

where N is the mean photon number of the reservoir which is 
related to the temperature of the substract, and C is the damp- 
ing constant. Even in our case where the Hamiltonian and 
Liouvilian are quadratic on quadrature operators, the solution 
of the problem is quite involved due to matrix exponentiation 

(14) 

and integration in (ITSJI. Next subsection will be devoted to 



the discussion of some particular regimes where the solution 
is feasible either analytically or numerically. 



IV. RESULTS 



C. Interaction with Environment 

The most general linear Markovian evolution of a density 
operator p preserving its positivity and the Gaussianity is de- 
scribed by the master equation (written in the Lindblad form): 



dp 
di 



H,p 
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with a quadratic Hamiltonian such as (|5]l and linear Lindblad 
superoperators Lk 



(16) 



where J is given by ^ and Afc is a complex vector By defin- 
ing the superoperators D = ReT, and F = J — ImT) 



where T = X^feli -^fc^L show that the covariance 

matrix (|9]l obeys now the following equation of motion ETl : 



(17) 



whose solution satisfying the initial condition 7(t = 0) = 7o 
is given bv Il27ll28l 



7(i) = * 70 e 



^^*+ rdi'er(*-*')Der^ (*-*'). (ig) 
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When thinking in terms of Wigner functions in phase space, 
all environment influence is represented by T, with its real 
part (contained in F) responsible for dissipation (the XH term 
alone is the symplectic evolution) while ImT is responsible 
for diffusion. 

In our system, the nanoscale motion of the NEMS is cer- 
tainly more susceptible to loss and decoherence than the mo- 
tion of the trapped ions. This happens because the former 
is subjected to coupling to phonons in the substrate where the 
nanobeam is fixed, while the ions are kept trapped by means of 
electromagnetic potentials in low pressure environments (al- 
most no residual gases) where collisions are minimized. Then, 
we will consider 



Li = A/C(iV + 1) a and L2 = 



(19) 



In what follows we present our results concerning entan- 
glement dynamics in the hybrid ion-NEMS-ion system pre- 
sented in last section. We will be considering physically mo- 
tivated cases consisting of closed and open system dynamics 
with zero or finite temperature T. In particular, we will be 
considering identical traps and ion species leading to equal 
frequencies for the ions vi = V2 — v and also equal ion- 
NEMS separation implying in ni = HI2 — n. In this symmet- 
ric scenario, the set of eigenvalues (or the spectrum) of XH is 
Specc(J'H) = {?7±,7r±,p±}, with 



7?± 

7r± 
P± 



±i 
±i 

71 



^2 + + 



a;2 — w^l + Sn?Luuj- 



(21) 



where we defined — v ^VL. 

It is the K-dependence of the quartet {7r±, p±\ which will 
determine the type of dynamics. Complex conjugate pairs 
of eigenvalues in SpeCj;^(J'H) indicate rotations whereas a 
real pair of eigenvalues with opposite signs indicates hyper- 
bolic movement in the phase space, the same as squeezing. 



From (21 1 it is easy to see that the first two eigenvalues r\± 
are always purely imaginary. Besides, if k < \Juji{y — fJ), 
Spcc£(JH) will consist only of pure imaginary eigenvalues 
indicating rotation in phase space. On the other hand, if 
K > \J lli{v — il), we obtain a "mixed" dynamics of rotations 
and "squeezing" which is then capable of generating some 
entanglement in the system. The same dynamical behavior 
can be obtained with more compact equations if we establish 
bj = V — Since v >> 51, in general, this corresponds to 
the limit of approximate resonance of natural frequencies of 
the ion vibrational motion and NEMS. In this case the above 
spectrum simplifies to 

Specc(J'H) = |±iw+, ±i-\/w(a; + k), ±i^a;(a; — k)| , 

(22) 

such that by taking k > a;, one will get a mixed dynamics 
(rotationsH-squeezing). We recall that generating squeezing is 
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a necessary condition to have entanglement from a separa- 
ble initial coherent or thermal state following Gaussian evo- 
lutions. Therefore, we shall now fix k > uj. The matrix ele- 



ments of Et necessary to obtain ( 12 1 for this spectrum can be 
found in appendix [A] 



A. Closed system with T = 

In this subsection, we will be considering an initial pure 
state consisting of the product of coherent states for the ions 
and the NEMS, this leads to the initial CM 70 = ^le, where 
Ig is the 6x6 identity matrix. We then evolve this initial co- 
variance matrix according to 



12 1. In order to facilitate the 



interpretation of the matrix elements, we arrange it in sectors 
representing each system and it is given by 



lit) 
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(23) 



where 




(24) 



The submatrices are shown in Appendix[A| It is worthwhile to 
notice the independence of the system CM ( 23 1 on cj+ which 



is the only parameter in ( 22 1 that contains information on the 
ion-ion coupling constant 51. This independence may be un- 
derstood as a consequence of the fact that uj± appears solely in 
the imaginary eigenvalues (rotations) not affecting the second 
moments. We are then in a regime where all entanglement is 
due to the coupling of the ions to the NEMS and not due the 
coupling between the ions. 

In figure |2j the NEMS -induced entanglement between the 
ions is presented as a function of time. This plot shows that the 
stronger the NEMS-ion coupling k, the stronger the entangle- 
ment between the ions. This coupling constant is externally 
controlled by the experimentalist and cj is fixed by the NEMS 
fabrication. Basically, according to Q, n can be enhanced by 
mere increasing of the gate voltage on the NEMS. It should 
be stressed that there is a lower bound on n for this creation of 
entanglement to happen (k must be stronger than w). Conse- 
quently, figure |2]illustrates a central result in this paper Two 
trapped ions initially uncorrected and prepared in coherent 
states can become entangled by interacting with a nanoelec- 
tromechanical resonator (also prepared in a coherent state) as 
soon as the ion-NEMS coupling achieve a certain value, and 
this can be controlled by external voltage gate on the NEMS 
device. 



B. Closed system with T 7^ 

A thermal state is a thermal equilibrium state of a quantum 
oscillator with frequency u) in contact with a thermal reservoir 




Figure 2: Logarithmic negativity for the ions when the system is 
initially prepared in a product of coherent states. Two coupling con- 
stants between the ions and the NEMS were considered: k = 1 
(solid) and k = 3 (dashed). Also, uj — 0.5 has been considered. 



at temperature T (canonical ensemble). Its density matrix be- 
comes diagonal in the Fock state basis |m): 



1 



1 ^ \n 

771 — 
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where the mean phonon number is 



|to)(to|, 



> 



(25) 



(26) 



and /j, is the Boltzman constant. The CM of such state is 
written as 'Yt •= + 5)12- By defining a — 2n+l and still 
considering the ions in a product of coherent states, the initial 
CM of the total system is 7q = iDiag(Q;l2, 12, 12) and its 
evolution, following ( [T2| and ( [22] i, will be 



7a (0 = 




with 



7al 



A1 Cj" 



(27) 



Cf^ Af 



(28) 



It has the same structure of equation ( 23 1 but now Cf and Af 
are lengthy and cumbersome expressions which we will not 
show explicitly. In the limit of T — > (a — > 1), we have 
checked analytically that 7a(t) — > 7(i) which is given by 

Let us now investigate the entangling power of a thermal 
NEMS. Firstly, let us study the ion-ion dynamics for two dif- 
ferent temperatures as shown in figure |3] We can see that 
by increasing the temperature (increasing a), two interesting 
trends appear. On one side, we can clearly see that the en- 
tanglement generated between the ions for a finite tempera- 
ture is smaller than in the T — case. In a certain way, this 
is expected since thermal fluctuations in the NEMS tend to 
turn it into a classical object. On the other hand, we can see 
that for a finite temperature there are collapses and revivals 
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Figure 3: Logarithmic negativity for the ions when the NEMS is pre- 
pared in a thermal state characterized by a and the ions are initially 
prepared in a product of coherent states. Thin lines correspond to 
zero temperature (a — 1) and thick lines correspond to a finite tem- 
perature (a — 5). For the NEMS-ions coupling we used ft — 1 
(solid) and k — 3 (dashed). Also, cj — 0.5 has been considered. 



of entanglement in the dynamics. Also, we can see that an 
increase in temperature leads to a decrease in the maximal 
amount of entanglement generated between the ions. Numeri- 
cal exploitation shows that there is not a temperature threshold 
above which no ion-ion bipartite entanglement is generated 
in the regime considered here. However, for higher tempera- 
tures the negativity remains non-zero (but very small) only for 
very short time intervals due to sudden death of entanglement 
and revivals. These trends can be seen in Fig|4] In general, 
the NEMS' entangling power degrades asymptotically as a 
grows, i.e., as the NEMS becomes more classical. 



vestigate the resilience of bi- and tripartite entanglements in 
the system. To treat the genuine tripartite entanglement in our 
system, we use the methods explained in Section [TlIB| In all 
the following plots, solid lines refer to the bipartite ion-ion 
entanglement, dashed lines refer to the bipartite NEMS-ion 
entanglement and dotted lines refer to the tripartite entangle- 
ment. Also, the initial system preparation is a product of co- 



herent states for each system, just like in Section IV A We 
would also like to emphasize that the bipartite entanglement 
between each ion and the NEMS are identical in the frequency 
regimes considered here, as previously explained. 

In figure |5] we present the open system time evolution of 
bi- and tripartite entanglements for two different NEMS-ion 
couplings K. The general behavior does not change with k, 
except the maximal value of the entanglements which clearly 
increases with k. From this plot, it is quite clear that the tri- 
paitite entanglement is much more resilient than the bipartite 
ones to losses in the NEMS. The latter goes to zero much later 
than the former as we can check numerically. 




Figure 4: Logarithmic negativity for the ions as a function of time 
t and a which is proportional to the temperature of the NEMS (see 
main text). The ions were initially prepared in a product of coherent 
states. The other parameters are uj — 0.5 and k, — 3. 



C. 



Open System Dynamics and Tripartite Entanglement 
Resilience 



We now study the case presented in Section lTlI C| where we 
briefly reviewed how dissipation in the NEMS can be taken 
into account in this Gaussian system. Our main goal is to in- 



Figure 5: Open system dynamics for bi- and tripartite entanglements 
in the system consisting of one NEMS coupled to two trapped ions. 
Thin lines refer to k = 1 and thick lines to k = 1.5. The other 
system parameters are = 0.01 and uj — 0.5. The thermal number 
of photons in the reservoir is iV — 4.5. Dotted lines corresponds to 
T012, dashed lines to Noi and continuous lines to A^i2. 

In order to further understand this apparent resilience of the 
tripartite entanglement when compared to the bipartite ones, 
we now present their time evolution for two different NEMS 
decay rates, keeping the reservoir temperature constant. This 
is presented in figure [6] As could be expected, the stronger 
the decay rates the stronger the suppression of entanglement. 
However, it is still remarkable that the tripartite entanglement 
takes a much longer time than the bipartite ones to go to zero, 
signalizing again a pronounced robustness to coupling to the 
environment. 

Now we analyze the thermal effect on the open system dy- 
namics. The results for two different temperatures is shown 
in figure [7] Although thermal noise severely degrades entan- 
glement, as expected, it is still noticeable that tripartite entan- 
glement is more robust then the bipartite ones. An increase in 
temperature is more destructive to the bipartite entanglement 
in our system than to the tripartite one in the sense that the lat- 
ter goes to zero much sooner than the former. We can see, for 
instance, that for iV = 15 no ion-ion bipartite entanglement 
appears, this shows that thermal fluctuations can erase this ef- 
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Figure 6: Open system dynamics for bi- and tripartite entanglements 
in tiie system for two different NEMS's decay constants. Thin lines 
refer to ( — 0.01 and thick lines to ( — 0.02. The other system 
parameters are n — 1 and u — 0.5. The thermal number of photons 
in the reservoir is iV = 4.5. Idem fig|5] 



feet of NEMS-indueed ion-ion entanglement. Actually there 
is a threshold for the temperature of the bath, numerically we 
found that above N = 14.2, no ion-ion is found (with preci- 
sion 10^^). 



B 1.0 



















//^'»*. "'^^'^Vr 

#"N / \ ,*•».'..— 


....:\ 







t 



Figure 7: Open system dynamics for bi- and tripartite entanglements 
in the system for two different temperatures. Thin lines refer to TV = 
4.5 and thick lines to iV = 15. The other system parameters are 
fc = 1, cj = 0.5 and ( = 0.01. Idem fig|5] 



D. Conclusions 

We have presented some results concerning entanglement 
generation in a tripartite system consisting of two atomic ions 
and a nanoscale mechanical resonator We showed that the lat- 
ter is capable of inducing entanglement in the former and we 
studied the features of such generation under ideal and non- 
ideal conditions. In the ideal case we show analytically the 
precise regime where the system becomes entangled or not. 
The amount of entanglement is controlled by the gate volt- 
age applied on NEMS. Especially, we found that there is not 
a temperature threshold above which no bipartite ion-ion en- 
tanglement is created when we are dealing with the ideal case 
(no dissipation), there is only a asymptotically degradation of 
entanglement as the temperature of the NEMS' initial state 
grows, surprisingly a distinct effect occurs when the NEMS 
is coupled to an environment: correlations between them take 



place and actually a temperature threshold appears. We have 
also found that the tripartite entanglement is much more re- 
silient to losses in the NEMS than the bipartite ones as the 
latter goes to zero much later than the former. We believe 
our studies may be useful for deepen our understanding about 
the interaction of atomic and nanoscale systems, especially 
in what concerns the appearing and destruction of nonlocal 
quantum correlations under ideal and nonideal conditions. 
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Appendix A: Some analytical matrix elements 

In this appendix, we present the explicit analytical form of 
some matrix elements needed to obtain results in the main 
body of this manuscript. 



The matrix elements of Ef ( 13 i relative to the spectrum ( 22 1 
are given by 



Eni = 


Et44 


= ^(coshwi + 


COSLUt), 


Etl2 = 


Et21 


— Et45 — E454 


~ ^ (cosh ujt — COS ult) 


Etl3 = 


Et31 


= Et4g — Eig4 


— ~^tl2, 


Et22 = 


Et55 


= j(coshwt + 


COS Ult + 2 COS Ld+t) , 


Et23 — 


Et32 


— Et56 = Etg5 





~ |(2cosw+t — coswi — coshwt), 



Et33 
Etl4 

Etl5 
Etl6 
Et25 
Et26 
Et41 

Et42 
Et43 
Et52 
Et53 



Et66 — j(2cosw+t + cosujt + cosh (it), 



](^ sin Cut + " sinhojt), 

Et24 = 



- ^(-^smujt + ^ sinhwt), 

4 V LJ 



— Et34 — — Ef]^5, 



4^u: 

1 / LJ 

1 



. ( " sill Ult + sinh ut — 2 sin ui+t) , 



Et3g = sindit + 1^ sinhwt + 2sinw+i) 

t35 

^{-^sinuit + sinhtji), 

2 ^ OJ LJ 

Et5i = ^(5sintDi+gsinh£:;i), 

Et61 — ~Ef42, 

Et63 = ^sinujt- 2 sin uj+t), 



4 ^ 

Et62 = i(S '^i^'^* ~ S sinhwt - 2sinw+t), 
where a) = \/uj{k — a;), Co = y^uj{K + uj) and uj^ := fl + v. 
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The submatrices necessary to write the covariance matrix 
([23]l are given by 



7^ = 



1 / 1 + a(t) + bit) c{t) + d{t) \ _ 

2 1^ c{t) + d{t) l + a'{t) + h'{t) j ' 

^ ^ ( l + \[a{t) + h{t)] \[c{t) + d{t)] 
' \ l[c{t) + d{t)] l + \[a'{t) + h'{t)] 

C ^ AI a{t) + b{t) c{t)+d{t) \ . 
' 2 \ c{i) + d{t) a'{t) + b'[t) ) ' 

^ a{t)~b{t) c{t)~d{t) \ 

4 I c{t) - d{t) a'(t) - b'{t) j ' 



where 



a{t) = sin^ut, (Al) 

^2 _ ^2 

a'{t) = r— sin^wi, (A2) 

2a; ^ 

LO^ I ^2 

b{t) = — — -— sinh^ wi, (A3) 
2a; ^ 



= ^ sinli^tl;^, (A4) 



c(t) = , _ sin2(Ii^, (A5) 
4a; a; 

o;*^ I ^2 

= sinh2aii. (A6) 
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